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Abstract 

We give an explicit recursive formula for the all £-loop integrand for scattering amplitudes 
in Af = 4 SYM in the planar limit, manifesting the full Yangian symmetry of the theory. This 
generalizes the BCFW recursion relation for tree amplitudes to all loop orders, and extends 
the Grassmannian duality for leading singularities to the full amplitude. It also provides a new 
physical picture for the meaning of loops, associated with canonical operations for removing 
particles in a Yangian-invariant way. Loop amplitudes arise from the "entangled" removal of 
pairs of particles, and are naturally presented as an integral over lines in momentum-twistor 
space. As expected from manifest Yangian invariance, the integrand is given as a sum over 
non-local terms, rather than the familiar decomposition in terms of local scalar integrals with 
rational coefficients. Knowing the integrands explicitly, it is straightforward to express them in 
local forms if desired; this turns out to be done most naturally using a novel basis of chiral, tensor 
integrals written in momentum-twistor space, each of which has unit leading singularities. As 
simple illustrative examples, we present a number of new multi-loop results written in local form, 
including the 6- and 7-point 2-loop NMHV amplitudes. Very concise expressions are presented 
for all 2-loop MHV amplitudes, as well as the 5-point 3-loop MHV amplitude. The structure 
of the loop integrand strongly suggests that the integrals yielding the physical amplitudes are 
"simple", and determined by IR-anomalies. We briefly comment on extending these ideas to 
more general planar theories. 
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1 The Loop Integrand for J\f = 4 SYM Amplitudes 



Scattering amplitudes in gauge theories have extraordinary properties that are completely invisible 
in the textbook formulation of local quantum field theory. The earliest hint of this hidden structure 
was the remarkable simplicity of the Parke- Taylor formula for tree-level MHV amplitudes [lj[2]- 
Witten's 2003 proposal of twistor string theory [3] gave a strong impetus to rapid developments in 
the field, inspiring the development of powerful new tools for computing tree amplitudes, including 
CSW diagra ms [I] and BCFW recursion relations [5jj8] . At one-loop, very efficient on-shell methods 
now exist [9||10| and at higher- loop level generalizations of the unitarity-based method 11-14 



have made a five-loop computation possible 15 , which should soon determine the five-loop cusp 
anomalous dimension [16] . 

The BCFW recursion relations in particular presented extremely compact expressions for tree 
amplitudes using building blocks with both local and non-local poles. In a parallel development, an 
amazing hidden symmetry of planar N = 4 SYM — dual conformal invariance — was noticed first in 



multi-loop perturbative calculations [17| and then at strong coupling 18 , leading to a remarkable 
connection between null-polygonal Wilson loops and scattering amplitudes 18-27 . It was quickly 
realized that the BCFW form of the tree amplitudes manifested both full superconformal and dual 
superconformal invariance, which together close into an infinite-dimensional Yangian symmetry 



algebra 28 . Understanding the role of this remarkable integrable structure in the full quantum 
theory, however, was clouded by the IR-divergences that appear to almost completely destroy the 
symmetry at loop-level, leaving only the anomalous action of the (Bosonic) dual conformal invariance 



29-32 



1.1 Grassmannian Duality for Leading Singularities 



In 33 , a strategy for making progress on these questions was suggested. The idea was to find 
objects closely associated with scattering amplitudes which are completely free of IR-divergences; 
the action of the symmetries would be expected to be manifest on such objects, and they would 
provide "data" that might be the output of a putative dual theory of the S-Matrix. 

The leading singularities of scattering amplitudes are precisely objects of this sort. Thinking of 
loop amplitudes as multi-dimensional complex integrals, leading singularities arise from performing 
the integration not on the usual non-compact 'contours' over all real loop-momenta, but on compact 
contours 'encircling' isolated (and generally complex) poles in momentum space. As such, they 
are free of IR-divergences and well-defined at any loop order, yielding algebraic functions of the 
external momenta. Leading singularities were known to have strange inter-relationships and satisfy 
mysterious identities not evident in their field-theoretic definition. Morally they are also expected 
to be Yangian- invariant, although even this is not completely manifeslQ Thus leading singularities 
seem to be prime candidates for objects to be understood and computed by a dual theory. 



1 Indeed we will give a proof of this basic fact in the next section; a different argument for the same result is given 
in [Ml. 
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Such a duality was proposed in [33] , connecting leading singularities of color-stripped, n-particle 
N fc MHV scattering amplitudes in Af = 4 SYM to a simple contour integral over the Grassmannian 

y ^ k{Z) = vol[GL(fc)] J (l---k)(2---k+l) aa --(n---k-l) 2 s4[4 (CaaZ a ). (1) 

Here a = 1, • • • , n labels the external particles, and Z a are variables in CP 3 ' 4 . The original formu- 
lation of this object worked with twistor variables W a = (W a \rj a ), and was given as £ nj fc + 2(W) = 
3^n,fc+2(W). This was quickly realized [35] to be completely equivalent to a second form in momen- 
tum twistor space [36], with £ ni fc + 2(A, A, rj) = M-J^ffy x y n ,k{Z). Here the variables Z a = (Z a \rj a ) 
are the "momentum- twistors" introduced by Hodges [37] , which are the most natural variables with 
which to discuss dual superconformal invariance. Furthermore, these momentum twistors are sim- 
ple algebraic functions of the external momenta, upon which scattering amplitudes conventionally 
depend^ 

Since the Grassmannian integral is invariant under both ordinary and dual superconformal 
transformations, it enjoys the full Yangian symmetry of the theory, as has been proven more directly 
in 38 . In fact, it has been argued that these contour integrals in G(k,n) generates all Yangian 
invariants, P 



41 
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Leading singularities are associated with residues of the Grassmannian integral. Residue the 
orems 



43 imply many non-trivial and otherwise mysterious linear relations between leading sin- 
gularities. These relations are associated with important physical properties such as locality and 
unitarity |33|. 

Further investigations [44] identified a new principle, the Grassmannian "particle interpreta- 
tion", which determines the correct contour of integration yielding the BCFW form of tree ampli- 
tudes |45| . Quite remarkably, a deformation of the integrand connects this formulation to twistor 
string theory [44 46 ,47]. Furthermore, another contour deformation produces the CSW expansion 



of tree amplitudes 48 , making the emergence of local space-time a derived consequence from the 
more primitive Grassmannian starting point. 

The Grassmannian integral and Yangian-invariance go hand-in-hand and are essentially syn- 
onymous; indeed, the Grassmannian integral is the most concrete way of thinking about Yangian 
invariants, since not only the symmetries but also the non-trivial relationship between different 
invariants are made manifest; even connections to non-manifestly Yangian-invariant but important 
physical objects (such as CSW terms) are made transparent. 

Given these developments, we are encouraged to ask again: is there an analogous structure 
underlying not just the leading singularities but the full loop amplitudes? Does Yangian-invariance 
play a role? And if so, how can we see this through the thicket of IR-divergences that appear to 
remove almost all traces of these remarkable symmetries in the final amplitudes? 



To quickly establish notation and conventions, the momentum of particle a is given by p„ = x^ +1 — iJJ, and 
the point x£ in the dual co-ordinate space is associated with the line (Z a -i Z a ) in the corresponding momentum- 
twistor space. This designation ensures that the lines (Z a -i Z a ) and (Z a Z a +i) intersect, so that correspondingly, 
x^ +1 — x% = p a is null. (Bosonic) dual-conformal invariants are made with 4-brackets {abed} — euklZ^Z^ Zfr Z^ . 
An important special case is (i— 1 i j—1 j) — (i— 1 j)(xj—Xi) 2 ; 2-brackets (ij) are computed using the upper-two 

components of Zi, Zj and cancel out in dual-conformal expressions. For more detail see 

for 



The residues of G(k, n) are Yangian-invariant for generic momenta away from collmear limits. See 
important discussions of the fate of Yangian invariance in the presence of collinear singularities. 
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1.2 The Planar Integrand 

Clearly, we need to focus again on finding well-defined objects associated with loop amplitudes. 
Fortunately, in planar theories, there is an extremely natural candidate: the loop integrand itself! 

Now, in a general theory, the loop integrand is not obviously a well-defined object. Consider the 
case of 1-loop diagrams. Most trivially, in summing over Feynman diagrams, there is no canonical 
way of combining different f-loop diagrams under the same integral sign, since there is no natural 
origin for the loop-momentum space. The situation is different in planar theories, however, and 



this ambiguity is absent. This is easy to see in the dual x-space co-ordinates 17 . The ambiguity 
in shifting the origin of loop momenta is nothing other than translations in x-space; but fixing the 
x\, . . . ,x n of the external particles allows us to canonically combine all the diagrams. Alternatively, 
in a planar theory it is possible to unambiguously define the loop momentum common to all diagrams 
to be the one which flows from particle "1" to particle "2". 

At two-loops and above, we have a number of loop integration variables in the dual space 
x, y, . . . , z, and the well-defined loop integrand is completely symmetrized in these variables. 

So the loop integrand is well-defined in the planar limit, and any dual theory should be able to 
compute it. All the symmetries of the theory should be manifest at the level of the integrand, only 
broken by IR-divergences in actually carrying out the integration — the symmetries of the theory 
are broken only by the choice of integration contour. 

1.3 Recursion Relations for All Loop Amplitudes 

Given that the integrand is a well-defined, rational function of the loop variables and the external 
momenta, we should be able to determine it using BCFW recursion relations in the familiar way[^J 
At loop-level the poles have residues with different physical meaning. The first kind is the analog of 
tree-level poles and correspond to factorization channels. The second kind has no tree-level analog; 
these are single cuts whose residues are forward limits of lower- loop amplitudes. Forward limits are 
naively ill-defined operations but quite nicely they exist in any supersymmetric gauge theory, as 
was shown to one-loop level in [50] . There it was also argued that forward limits are well-defined to 
higher orders in perturbation theory in J\f = 4 SYM. In principle, this is all we need for computing 
the loop integrand in M = 4 SYM to all orders in perturbation theory. However, our goal requires 
more than that. We would like to show that the integrand of the theory can be written in a form 
which makes all symmetries — the full Yangian — manifest. The Yangian-invariance of BCFW terms 
at tree-level becomes obvious once they are identified with residues of the Grassmannian integral, 
we would like to achieve the same at loop-level. 

This is exactly what we will do in this paper. We will give an explicit recursive construction 
of the all-loop integrand, in exact analogy to the BCFW recursion relations for tree amplitudes, 
making the full Yangian symmetry of the theory manifest. 

The formulation also provides a new physical understanding of the meaning of loops, associated 
with simple operations for "removing" particles in a Yangian-invariant way. Loop amplitudes are 
associated with removing pairs of particles in an "entangled" way. We describe all these operations in 



We note that 49 have conjectured that the loop amplitudes can be determined by CSW rules, manifesting the 



superconformal invariance of the theory. 
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momentum-twistor space, since this directly corresponds to familiar momentum-space loop integrals; 
presumably an ordinary twistor space description should also be possible. 

As is familiar from the BCFW recursion relations at tree- level, the integrand is expressed as a 
sum over non-local terms, in a form very different than the familiar "rational function x scalar inte- 
gral" presentation that is common in the literature. Nonetheless, the Yangian-invariance guarantees 
that every term in the loop amplitude has Grassmannian residues as its leading singularities. 

The integrands can of course be expressed in a manifestly-local form if desired, and are most 



naturally written in momentum-twistor space [51,52 . As we will see, the most natural basis of local 



integrands in which to express the answer is not composed of the familiar scalar loop-integrals, but 
is instead made up of chiral tensor integrals with unit leading-singularities, which makes the physics 
and underlying structure much more transparent. 

Of course the integrand is a well-defined rational function which is computed in four-dimensions 
without any regulators. The regularization needed to carry out the integrations is a very physical 
one, given by moving out on the Coulomb branch [53] of the theory. This can be beautifully 
implemented, both conceptually and in practice, with the momentum-twistor space representation 



of the integrand 51 , 52 



Quite apart from the conceptual advantages of this way of thinking about loops, our new for- 
mulation is also completely systematic and practical, taking the "art" out of the computation of 
multi-loop amplitudes in M = 4 SYM. As simple applications of the general recursive formula, 
we present a number of new multi-loop results, including the two- loop NMHV 6- and 7-particle 
integrands. We also include very concise, local expressions for all 2-loop MHV integrands and for 
the 5-particle MHV integrand at 3-loops. All multiplicity results for the so-called "parity even" 
part of two-loop amplitudes in the MHV sector were obtained by Vergu in (54), extending previous 



work done for 5-particles 15 and 6-particles 25 , 55] in dimensional regularization. The "parity 
even" part of the 6-particle amplitude in dimensional regularization has been computed in work in 
progress by Kosower, Roiban, and Vergu [56]. Complete integrands have been computed at two- 



loop order for 5-particles in [15] using D-dimensional unitarity and for 5- and 6-particles in [55,57 
using the leading singularity technique developed in [57, 58 . Also using the leading singularity 



technique, the 5-point 3-loop integrand was presented in [59]. Combining D-dimensional unitarity 
with a generalization of quadruple cuts to higher loop order [58] , a method called maximal cuts 
was introduced in 15 and used for the computation of the 4-point 5-loop integrand. The 4-point 



amplitude integrand at I = 2, 3, 4 loop- level were computed in [60], [61], and 62 , respectively. The 
method to be used in this paper is, however, very different both in philosophy and in practice from 
the leading singularity or generalized unitarity approaches. 

In this paper, we give a brief and quite telegraphic outline of our arguments and results; we will 
present a much more detailed account of our methods and further elaborate on many of the themes 



presented here in upcoming work 63 . In section [2J we describe a number of canonical operations 
on Yangian invariants — adding and removing particles, fusing invariants — that generate a variety 
of important physical objects in our story. In section [3] we describe the origin of Yangian-invariant 
loop integrals as arising from the "hidden entanglement" of pairs of removed particles. In section [4] 
we describe the main result of the paper: a generalization of the BCFW recursion relation to all loop 
amplitudes in the theory, and discuss some of its salient features through simple 1-loop examples. 
In section [5] we set the stage for presenting loop amplitudes in a manifestly local form by describing 
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the most natural way of doing this in momentum- twistor space. In section [6] we present a number 
of new multi-loop integrands computed using the recursion relation and translated into local form 
for the convenience of comparing with known results where they are available. We conclude in 
section [7] with a discussion of a number of directions for future work. We discuss indications that 
not only the integrands but also the loop integrals should be "simple" . The idea of determining the 
loop integrand for planar amplitudes is a general one that can generalize well beyond maximally 
supersymmetric theories with Yangian symmetry, and we also very briefly discuss these prospects. 



2 Canonical Operations on Yangian Invariants 

As a first step towards the construction of the all-loop integrand for M = 4 SYM in manifestly 
Yangian form, we study simple operations that can map Yangian invariants Y n k{Z\, ■ ■ ■ ,2 n ) to 
other Yangian invariants. In this discussion it will not matter whether the Z's represent variables 
in twistor-space or momentum-twistor space; we will simply be describing mathematical operations 
that mapping between invariants. Combining these operations in various ways yields many objects 
of physical significance |63| . The same physical object will arise from different combinations of 
these operations in twistor-space vs. momentum-twistor space; we will content ourselves here by 
presenting mostly the momentum-twistor space representations. 

As mentioned in the introduction, understanding these operations is not strictly necessary if we 
simply aim to find a formula for the integrand. The reason is that the BCFW recursion relations we 
introduce in section|4]can be developed independently for theories with less supersymmetry, which do 
not enjoy a Yangian symmetry. Our insistence in keeping the Yangian manifest however will pay off 
in two ways. The first is conceptual: the Yangian-invariant formulation will introduce a new physical 
picture for meaning of loops. The second is computational: the Yangian-invariant formulation gives 
a powerful way to compute the novel "forward-limit" terms in the BCFW recursions in momentum- 
twistor space, using the Grassmannian language. 

We will begin by discussing how to add and remove particles in a Yangian-invariant way. One 
motivation is an unusual feature of the Grassmannian integral-the space of integration depends on 
the number re of particles. It is natural to try and connect different re's by choosing a contour of 
integration that allows a "particle interpretation" , by which we mean simply that the variety defining 
the contour for the scattering amplitudes of (re + 1) particles differs from the one for n particles only 
by specifying the extra constraints associated with the new particle [441. Following this "add one 



particle at a time" -guideline completely specifies the contour for all tree amplitudes 44,47], along 
the way exposing a remarkable connection with twistor string theory [3 46 64 66] . As we will see in 
this paper, loops are associated with interesting "entangled" ways of removing particles from higher- 
point amplitudes. We will then move on to discuss how to "fuse" two invariants together. Using 
these operations we demonstrate the Yangian invariance of all leading singularities, and discuss the 
important special case of the "BCFW bridge" in some detail. 
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2.1 Adding Particles 

Let us start with a general Yangian-invariant object 

Y n> k(2i, . . . , Z n ). (2) 

We will first describe operations that will add a particle to lower-point invariants to get higher- 
point invariants known as applying "inverse soft factors" [67], which are so named because taking 
the usual soft limit of the resulting object returns the original object. This can be done preserving 
k or increasing k \— > k + 1. We can discuss these in both twistor- and momentum- twistor space; for 
the purposes of this paper we will describe these inverse-soft factor operations in momentum-twistor 
space. 

The idea is that there are residues in G(k, n) which are trivially related to residues in G(k, n — 1) 
or G(k — l,n — 1). The ^-preserving operation Y n _\^ h> Y n ^ is particularly simple, being simply 
the identification 

Y n ^Z\i ■ ■ ■ ,Zn-l,Zn) = Yn-l,k{Z\, ■ ■ . , Z„,_i); (3) 

that is, where we have simply added particle n as a label (but have not altered the functional 
form of Y in any way); thanks to the momentum-twistor variables, momentum conservation is 
automatically preserved. The /c-increasing inverse soft factor is slightly more interesting. There is 
always a residue of G(k, n) which has a C-matrix of the form 

* * ••• * * 1 \ 

* * _ (4) 

; ••. •■. ••• ••. ■.*:) 

Here, the non-zero elements in the top row, * * * * 1 correspond to particles 1, 2, (n — 2), (n — 1), n, 
and we have generic non-zero entries in the lower (k — 1) x (n — 1) matrix. The corresponding 
residue is easily seen to be associated with 

Y^ k (. . . ,Z n -i,Z n ,Z\, ...) = [n-2 n-1 nl2] x F n _ 1)fc _i(. . .,Z n -i,Zx, . . .) (5) 

where 

\ abcdc \ = S°\ 4 (r) a (bcde) + cyclic) 

(abcd)(bcde)(cdea){deab)(eabc) 

is the basic 'NMHV'-like i?-invariant|^] and the Z n —\^\ are deformed momentum twistor variables. 
The Bosonic components of the deformed twistors have a very nice interpretation: Z\ is sim- 
ply the intersection of the line (12) with the plane (n— 2 n— 1 n), which we indicate by writing 
Z\ = (n—2 n—1 n) (~) (1 2); and Z n _\ is the intersection of the line (n— 2 n— 1) with the plane (1 2 re), 
written Z n -\ = (n—2 n—1) f] (n 1 2). Fully supersymmetrically, we have 

Z\ = (n-2 n-1 n) fl (1 2) = Z x (2 n-2 n-1 n) + Z 2 (n-2 n-1 n 1); 
Z n _i = (n-2 n-1) D in 1 2) = Z n _ 2 (n-1 n 1 2) + Z n _i(n 1 2 n-2). 



J When two sets of the twistors are consecutive, these "_R-invariants" are sometimes written R r - S ,t = [r s-1 s t-1 t]. 
These invariants were first introduced in [26] in dual super-coordinate space. 
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2.2 Removing Particles 

We can also remove particles to get lower-point Yangian invariants from higher-point ones. This 
turns out to be more interesting than the inverse-soft factor operation, though physically one might 
think it is even more straightforward. After all, we can remove a particle simply by taking its soft 
limit. However, while this is a well-defined operation on e.g. the full tree amplitude, it is not a 
well-defined operation on the individual residues (BCFW terms) in the tree amplitude. The reason 
is the presence of spurious poles: each term does not individually have the correct behavior in the 
soft limit. 

Nonetheless, there are completely canonical and simple operations for removing particles in a 
Yangian-invariant way. One reduces fc i — >■ A; — 1, the other preserves k. The fc-reducing operation 
removes particle n by integrating over its twistor co-ordinate 

Y n-l,k-l{Z\, • • • Zn-l) = J d 3 ^Z n Y n ${Z\, . . . , Z n -i,Z n ). (8) 

This gives a Yangian-invariant for any closed contour of integration — meaning that under the Yan- 
gian generators for particles 1, . . . , n — 1, this object transforms into a total derivative with respect 
to Z n . This statement can be trivially verified by directly examining the action of the level-zero 
and level-one Yangian generators on the integral. It is also very easy to verify directly from the 
Grassmannian integral. Note that depending on the contour that is chosen, a given higher-point 
invariant can in general map to several lower-point invariants. 

The fc-preserving operation "merges" particle n with one of its neighbors, n — 1 or 1. For 
example, 

• • -Zn-i) = y n ,k(Zi, ■ ■ . ,Z n _i,Z n i->- Z n -i). (9) 

The Yangian-invariance of this operation is slightly less obvious to see by simply manipulating 
Yangian generators, but it can be verified easily from the Grassmannian formula. 

We stress again that these operations are perfectly well-defined on any Yangian-invariant object, 
regardless of whether the standard soft-limits are well defined. Of course, they coincide with the 
soft limit when acting on e.g. the tree amplitude. 

2.3 Fusing Invariants 

Finally, we mention a completely trivial way of combining two Yangian invariants to produce a new 
invariant. Start with two invariants which are functions of a disjoint set of particles, which we can 
label Y\{Z\, . . . , Z m ) and l^C^m+i, ■ • • , Z n ). Then, it is easy to see that the simple product 

Y'{Z 1 , ...,Z n ) = Y 1 (Z 1 , ...,Z m )x Y 2 (Z m+1 , ...,Z n ) (10) 

is also Yangian-invariant. Only the vanishing under the level-one generators requires a small com- 
ment. Note that the cross terms vanish because the corresponding level-zero generators commute 
and therefore the level-one generators cleanly splits into the smaller level-one generators. 
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2.4 Leading Singularities are Yangian Invariant 

Combining these operations builds new Yangian invariants from old ones; all of which have nice 
physical interpretations. An immediate consequence is a simple proof that all leading singularities 
are Yangian invariant. For this subsection only, we work in ordinary twistor space. Then we take 
any four Yangian invariants for disjoint sets of particles and we make a new invariant by taking the 
product of all of them, Y^Wi, . . - , W m )Y 2 (W m+ i, Wi)Y 3 (Wi+i, W p )Y±(W p+h ...,>%). We 
then "merge" m and m + 1, I and l + l, p and p + 1, and q with 1. We then integrate over m, l,p, q. 
This precisely yields the twistor-space expression for a "1-loop" leading singularity topology |68p9| . 




In the figure, a thick black line denotes the merging of the two particles at the ends of the line, and 
integrating over the remaining variable. The generalization to all leading singularities is obvious; 
for instance, starting with the "1-loop" leading singularity we have already built, we can use the 
same merge and integrate operations to build "2-loop" leading singularity topologies such as that 
shown below. 








We conclude that all leading singularities are Yangian invariant. Given that all Yangian invariants 



are Grassmannian residues, this proves (in passing) the original conjecture in 33 that all leading 
singularities can be identified as residues of the Grassmannian integral. 
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2.5 The BCFW Bridge 



A particularly important way of putting together two Yangian invariants to make a third is the 
"BCFW bridge" (7 70 , 71 , associated with the familiar "two-mass hard" leading singularities drawn 



below in twistor space 70 - 73 




n-1 



Here, the open and dark circles respectively denote MHV and MHV three-particle amplitudes, 
respectively. We remark in passing that the inverse-soft factor operations mentioned above are 
special cases of the BCFW bridge where a given Yangian invariant is bridged with an MHV three- 
point vertex (for the /c-preserving case) or an MHV three-point vertex (for the fe-increasing case). 

We will find it useful to also see the bridge expressed as a composition of our basic operations 
in momentum-twistor space, as 




2 1 




71 — 1 




2 1 In n-1 



This is a pretty object since it uses all of our basic operations in an interesting way. In the figure, 
the solid arrows pointing inward indicate that particle- "1" is added as an fc-increasing inverse 
soft factor on Yl, and j+1 is added as a /c-increasing inverse soft factor on Y#. We are also 
using the merge operation to identify the repeated "1" and "j+1" labels across the bridge. The 
internal line, which we label as Zj, is integrated over. The contour of integration is chosen to 
encircle the (n— 1 n 1 /)-pole from the [n—1 nil j+l]-piece of the inverse-soft factor on Yl, and 
the (1 / j+1 j)- and (I j+1 j j— l)-poles from the [1 / j+1 j j— l]-piece of the inverse soft factor 
on Yr. The deformation on Z n induced by the inverse-soft factor adding particle- 1 on Yl is of the 
form 

Z n <-+ Z n = Z n + zZ n -i, where (Z n ZiZjZj + x) = 0. (11) 

This is the momentum-twistor space version of the BCFW deformation, which corresponds to 
deforming X n , Ai in momentum-space. We remind ourselves of this deformation by placing the little 
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arrow pointing from n i— > n — 1 in the figure for the bridge. The momentum-twistor space geometry 
associated with this object is 




71—1 



which precisely corresponds to the expected BCFW deformation and the corresponding factorization 
channel. 

We leave a detailed derivation of this picture to [63] , but in fact the momentum-twistor structure 
of the BCFW bridge can be easily understood. Note that Yl, Yr have fc-charge Jcl, kn, while Yl<8>Yr 
has fc-charge fcz, + &r + 1; given that the Zi decreases the fc-charge by 1, we must start with Yj, 
and Yr and get objects with fc-charge (ki + 1) and (/cr + 1) on the left and right. This can be 
canonically done by acting with ^-increasing inverse soft factors; the added particle on Yl must be 
adjacent to n in order to affect a deformation on Z n . Finally, the data associated with the "extra" 
particles introduced by the inverse soft factor must be removed in the only way possible, by using 
the merge operation. Explicitly, the final result for Yl <S> Yr is 

BCFW 

(Yl ® Y R )(l,...,n) = [n-ln lj j+1] x Y R (l, . . . ,j,l) x Y L {l,j + 1, . . . ,n-l,n) (12) 

V BCFW / 

Wlth n = (n-1 n)f](j j+1 1), and I = (j j+1) f] (n-1 n 1). (13) 

Starting with the tree amplitude M ni fc,tree|^J the BCFW deformation Z n h > Z n + zZ n -\ can be 
used to recursively construct tree amplitudes in the familiar way: by writing, 

/dz 
— M nAtrec (z), (14) 

it is clear that the desired amplitude M n> k^ Tee (z) is obtained by summing-over all the residues of 
the RHS except the pole at origin z = 0. Notice that there is a non-zero pole at infinity in this 
deformation: as z — > oo, Z n — > Z n ^\ projectively, and so the tree amplitude gets a contribution 
from M n (Z\, . . . , Z n -\, Z n ) — > M n -\{Z\, . . . , Z n -\) ^ The pole at z — > oo corresponds to the term 

6 We remind the reader that we are working in momentum-twistor space, so that what we are calling Mtroo here is 
obtained after stripping off the MHV tree-amplitude factor from the full amplitude in momentum space. 

7 Note that z — > oo here does not correspond to going to infinity in the familiar momentum-space version of 
BCFW. The pole at infinity in ordinary momentum space here corresponds to a pole involving the infinity twistor 
(Z n (z) I Z\) = 0. Of course we do not expect such a pole to arise in a dual-conformal invariant theory, not only 
at tree-level, but at all-loop order, as will be relevant to our subsequent discussion. A direct proof of this fact, not 



assuming dual conformal invariance, should follow from the "enhanced spin-lorentz symmetry" arguments of 71 
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in the usual momentum-space BCFW formula using an MHV three-point vertex bridged with JVf n _i, 
which simply acts as a ^-preserving inverse-soft factor The remaining physical poles are of the form 
(i i+1 j Under Z n \- y Z n + zZ n -i, we only access the poles where (Z n (z)ZiZjZj+i) — > 0, 

and the corresponding residues are computed by the BCFW bridge indicated above, with Yl,Yr 
being the lower-point tree amplitudes. 



3 Loops From Hidden Entanglement 

Let's imagine starting with some scattering amplitude or Grassmannian residue, and begin remov- 
ing particles. The operation that decreases k in particular demands a choice for the contour of 
integration. If we remove particle Za by integrating over it as J cP^Za, it is natural to choose a 
T 3 -contour of integration for the Bosonic cPZa integral and compute a simple residu^J 

We can then proceed to remove a subsequent particle either by merging, or performing further 
integrals J cP^Zb and so on. In this way we will simply proceed from higher-point Grassmannian 
residues to lower-point ones. In particular, if these operations are performed on a higher-point tree 
amplitude, we arrive at lower-point tree amplitudes, and don't encounter any new objects. 

But we can imagine a more interesting way of removing not just one but a pair of particles. 
Consider removing particle A and subsequently removing the adjacent particle B. Instead of first 
integrating-out A and then B on separate T 3 's, let's consider an "entangled" contour of integration, 
which we will discover to yield, instead of lower-point Grassmannian residue, a loop integral. 

Consider as a simple example removing two particles from the 6-particle N 2 MHV = MHV tree 
amplitude, M6 ) 4 r g=o(1234AB) . Performing the cP^ija, d ^ 4 r]B integrals is trivial, and this gives 

(1234) 3 

d ZB (2Mz A ) {Mz A z B ) (4z a z b 1) (z a z b 12) (^123) (15) 

where we have chosen to label the Bosonic momentum twistors with lower-case z'a for later conve- 
nience. As we have claimed, on any closed contour, these integrals should give a Yangian-invariant 
answer. Indeed, computing the zb integral by residue on any contour leaves us with 

d 3 z A -, n n r (16) 

and computing any of the simple residues of this remaining za integral gives 1, which is of course 
the only Yangian invariant for MHV amplitudes. 

We will now see that starting with exactly the same integrand but choosing a different contour of 
integration yields, instead of "1", the 4-particle 1-loop amplitude. Geometrically, the points za, zb 
determine a line in momentum-twistor space, which is interpreted as a point in the dual x-space, or 
equivalently, a loop-integral's four- momentum. We will first integrate over the positions of za, zb 
on the line (AB), and then integrate over all lines (AB). 



8 Residues of rational functions in m complex variables are computed by choosing m polynomial factors fi's from 
the denominator and integrating along a particular T m -contour, i.e. the product of m circles given as the solutions 



of I/; | — e with e C 1 and near a common zero of the fi's. See 43 for more details 



12 



This contour can be described explicitly by parametrizing z a ,b as 



ZA = ( XA ) z B =( Xb ) 



(17) 



where x will be the loop momentum. The measure is 



d 3 z A d 3 z B = {\ A d\ A ){\Bd\ B ){^A^B) 2 d A x. 



(18) 



The X A , A_b integrals will be treated as contour integrals on CP x CP 1 , while the x-integral will be 
over real points in the (dual) Minkowksi space. 

Using that (zazb j-1 j) = (AaAb)(j-1 j)(x—Xj) 2 our integral becomes 



The factor (^234) is linear in the projective variable A^ while the factor (123z.b) is linear in X B . 
This implies that there is a unique way to perform the A^ and A# integrals by contour integration, 
which gives us 



This is precisely the 1-loop MHV amplitude! 

We have thus seen that, removing a pair of particles with this "entangled" contour of integration, 
where we first integrate over the position of two points along the line joining them and then integrate 
over all lines, naturally produces objects that look like loop integrals. 

There is a nicer way of characterizing this "entangled" contour that is also more convenient for 
doing calculations, let us describe it in detail. Given za, Zb, a general GL(2)-transformation on the 
2-vector (za, zb) moves A, B along the line (AB). Thus, in integrating over d 3 ZAd 3 zs, we'd like to 
"do the GL(2)-part of the integral first" to leave us with an integral that only depends on the line 




(19) 




(20) 



(AB): 




We can do this explicitly by writing 




(21) 



then 



d 3 z A d 3 z B = (c A dc A )(c B dc B )(cAC B ) 2 



d 4 Z A d 4 Z B 
vol[GL(2)] 



(22) 
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and our integral becomes — this time writing it out fully: 

~ d 4 Z A d*Z B ^ (1234) 3 

_vol[GL(2)] 



/ 



(AB12)(AB34)(AB41) 



(c A dc A ) (c B dc B ) 
(c A c B )(c A ip A }{c B ip B ) 



where 



(A234) 
(B234) 



(A123) 
(5123) 



(23) 



(24) 



The c A ,c B integral is naturally performed on a contour 'encircling' c A = ip A ,c B = ipB, yielding 
! — 1 More generally, if "234" and "123" in the definitions of ip A ,ip B were to be re- 

(ip A ip B ) = (Axyz)(Babc)-(Aabc}{Bxyz) = (AB\(abc)C\(xyz)) 



{■4>a^b) (AB23>(1234>- 

placed by arbitrary "a&c" and u xyz" 
where (a&c) f] (xyz) is the line corresponding to the intersection of the planes (abc) and (xyz). We 
are then left with 



d A Z A d A Z B 
vol[GL(2)] 



(1234)' 



(AB12)(AB23)(AB34)(AB41)' ^ 

where the integration region is such that the line (AB) corresponds to a real point in the (dual) 
Minkowski space-time. We recognize this object as the 1-loop MHV amplitude, exactly as above. 

We can clearly perform this operation starting with any Yangian invariant object Y[Z A , Z B ,Z\, . . 
which we will graphically denote as: 

A 
/ -B 




\ 1 

2 

and write as 

J Y\...,Z n ,Z A ,Z B ,Z x ,..\ (26) 

GL(2) 

This object is formally Yangian- invariant, in the precise sense that the integrand will transform 
into a total derivative under the action of the Yangian generators for the external particles. Of 
course, such integrals may have IR-divergences along some contours of integration, which is how 
Yangian-invariance is broken in practice. 

The usual way of writing the loop amplitudes as "leading singularity x scalar integral" ensures 
that the leading singularities of the individual terms are Yangian-invariant, but the factorized 
form seems very un-natural, and there is no obvious action of the symmetry generators on the 
integrand. By contrast, the loop integrals we have defined, as we will see, will not take the artificial 
"residue x integral" form, but of course their leading singularities are automatically Grassmannian 
residues. The reason is that a leading singularity of the (yli3)-integral can be computed as a simple 
residue of the underlying d^z A d^z B integral, which is free of IR divergences and guaranteed to 
be Yangian-invariant. 
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4 Recursion Relations For Arbitrary Loop Amplitudes 



Having familiarized ourselves with the basic operations on Yangian invariants, we are ready to 
discuss the recursion relations for loops in the most transparent way. The loop integrand is a 
rational function of both the loop integration variables and the external kinematical variables. Just 
as the BCFW recursion relations allow us to compute a rational function from its poles under a 
simple deformation, the loop integrand can be determined in the same way. Consider the Z-loop 
integrand M n ^/, and consider again the (supersymmetric) momentum-twistor deformation 

Z n i-+ Z n + zZ n -x. (27) 

Then 

/dz 
— M nA£ (z) (28) 

and we sum over all the residues of the RHS away from the origin, all of which can be determined 
from lower-point/lower-loop amplitudes. This recursion relation can be derived in a large class of 
theories and is not directly tied to N = 4 SYM or Yangian-invariance. However our experience 
with building Yangian-invariant objects will help us to understand (and compute) the terms in 
the recursion relations in a transparent way, and also easily recognize them as manifestly Yangian- 
invariant objects. 

As in our discussion of the BCFW bridge at tree-level, the pole at infinity is simply the lower- 
point integrand with particle n removed. All the rest of the poles in z also have a simple interpre- 
tation: in general, all the poles arise either from (Z n (z) Z\ Zj ^j+i) — > or ((AB) q Z n {z) Z\) — > 0, 
where {AB) q denotes the line in momentum twistor space associated with the g th loop-variable. 
The first type of pole simply corresponds to factorization channels, and the corresponding residue 
is computed by the BCFW bridges between lower- loop/lower-point amplitudes: 



1 




where ul + ur = n + 2, fct + /cr = k — 1, II + £r = L Note that we treat all the poles 
(including the pole at infinity) on an equal footing by declaring the term with j = 2 to be given by 
the fc-preserving inverse soft-factor acting on lower-point amplitude. 

This is the most obvious generalization of the BCFW recursion relation from trees to loops, but 
it clearly can't be the whole story, since it would allow us to recursively reduce loop amplitudes to 
the 3-particle loop amplitude, which vanishes! Obviously, at loop-level, a "source" term is needed 
for the recursive formula. 
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4.1 Single-Cuts and the Forward-Limit 



This source term is clearly provided by the second set of poles, arising from ((AB) q Z n (z) Z\) — > 0. 
For simplicity of discussion let's first consider the 1-loop amplitude. This pole corresponds to cutting 
the loop momentum running between n and 1, and is therefore given by a tree-amplitude with two 
additional particles sandwiched-between n, 1, with momenta q, —q, summing-over the multiplet of 
states running around the loop. These single-cuts associated with "forward-limits" of lower-loop 
integrands are precisely the objects that make an appearance in the context of the Feynman tree 
theorem |50|. The geometry of the forward limit is shown below for both in the dual x-space and 
momentum-twistor space: 



■x 7 




Here, between particles 5 and 1, we have particles 6, 7 with momenta q^j—q^, where q^ = 

is a null vector. In momentum-twistor space, the null condition means that the line (76) intersects 

(15), while in the forward limit both Z% and Z^ approach the intersection point (76) f] (15). 

In a generic gauge theory, the forward limits of tree amplitudes suffer from collinear divergences 



and are not obviously well-defined. However remarkably, as pointed out in 50 , in supersymmetric 
theories the sum over the full multiplet makes these objects completely well-defined and equal to 
single-cuts! 

Indeed, we can go further and express this single-cut "forward limit" term in a manifestly 
Yangian-invariant way. It turns out to to be a beautiful object, combining the entangled removal 
of two particles with the "merge" operation: 



n— 1 




Here a particle (n + 1) is added adjacent to A, B as a fc-increasing inverse soft factor, then A, B are 
removed by entangled integration. The GL(2)-contour is chosen to encircle points where both points 
A,B on the line (AB) are located at the intersection of the line (AB) with the plane (n— 1 n 1). 
Note that there is no actual integral to be done here; the GL(2) integral is done on residues and is 
computed purely algebraically. Finally, the added particle (n + 1) is merged with 1. 
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As in our discussion of the BCFW bridge, this form can be easily understood by looking at the 
deformations induced by the "1" inverse soft factors; the associated momentum-twistor geometry 
turns out to be 

(AB\ 

>Z A 



n—1 



exactly as needed. The picture is the same for taking the single cut of any Yangian- invariant object. 

Note that we were able to identify the BCFW terms in a straightforward way since the residues 
of the poles of the integrand have obvious "factorization" and "cut" interpretations. This is another 
significant advantage of working with the integrand, since as is well known, the full loop amplitudes 
(after integration) have more complicated factorization properties [74] . This is due to the IR 
divergences which occur when the loop momenta becomes collinear to external particles, when 
the integration is performed. 



4.2 BCFW For All Loop Amplitudes 

Putting the pieces together, we can give the recursive definition for all loop integrands in planar 
M = 4 SYM as 



n — 1 




To be fully explicit, the recursion relation is 
k Al, ...,n)= M n _i w(l, . . . , n - 1) 



+ J [AB n-l n 1] x M n+2 ,k+i,e-i{l, ■ ■ ■ ,n AB , A, B). 

GL(2) 



(29) 



where nz, + riR = n + 2, ki, + Hr = k — 1, £l + £r = £ and the shifted momentum (super-)twistors 
that enter are 



% = ( n_1 n ) H U j+1 1), 
n AB = (n-l n)f](ABl), 



Ij = (i J+1) f1 (n-l n 1) 
A= (AB)f](n-l n 1). 



(30) 
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Beyond 1-loop, it is understood that this expression is to be fully-symmetrized with equal weight in 
all the loop-integration variables (AB)f, it is easy to see that this correctly captures the recursive 
combinatorics. Recall again that GL(2)-integral is done on simple residues and is thus computed 
purely algebraically; the contour is chosen so that the points A, B are sent to (AB) f] (n—1 n 1). As 
we will show in [63], recursively using the BCFW form for the lower-loop amplitudes appearing in 
the forward limit allows us to carry out the GL(2)-integral completely explicitly, but the form we 
have given here will suffice for this paper. 



4.3 Simple Examples 



In [63], we will describe the loop-level BCFW computations in detail; here we will just highlight 
some of the results for some simple cases, to illustrate some of the important properties of the 
recursion and the amplitudes that result. We start by giving the BCFW formula for all one-loop 
MHV amplitudes. 

In this case the second line in the above formula vanishes, and the recursion relation trivially 
reduces to a single sum. To compute the NMHV tree amplitudes which enters through the third 
line, it is convenient to use the tree BCFW deformation Zb = Zb + zZa which leads to 



1— loop 
MHV 



d^Z A d^Z B 

vol[GL(2)] 



J Y^[ABjj+l 1] x (XJ[AB 1 

r i9\ 3 \ i<3 



i + ... 



(31) 



GL(2) 



where the omitted terms are independent of Zb and vanish upon Fermionic-integration. The GL(2) 
and Fermion integrals are readily evaluated, as explained above, reducing this to 



1— loop 
MHV 



d A Z A d A Z B 
vol[GL(2)] 



Kj 



(AB\(lii+l)f)(lj j+l))' 



{AB 1 i)(AB i i+l)(AB i+1 1)(AB 1 j)(AB j j+l)(AB j+1 1) 



.(32) 



This is the full one-loop integrand for MHV amplitudes. 

As expected on general grounds from Yangian-invariance, and also as familiar from BCFW re- 
cursion at tree-level, the individual terms in this formula contain both local and non-local poles. 
We will graphically denote a factor (ABxy) in the denominator by drawing a line (xy); the nu- 
merators of tensor integrals (required by dual conformal invariance) will be denoted by wavy- and 
dashed-lines — the precise meaning of which will be explained shortly. In this notation, this result is 



n—1 




i+1 




Notice that all the terms have 6 factors in the denominator, and hence by dual conformal invariance 
we must have two factors containing AB in the numerators. These are denoted by the wavy lines: the 
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numerator is {AB\(1 i i+1) 0(1 j j+l)) 2 = {(A 1 i i+l)(B 1 j j+l) - (B 1 i i+l)(A 1 j j+l)) 2 , 
where the power of 2 has been indicated by the line's multiplicity. 

Notice that when i + 1 = j, the numerator cancels the two factors (AB 1 j) 2 in the denominator: 
by a simple use of the Schouten identity it is easy to see that 

[(A 1 j-1 j)(B 1 j j+l) - {A 1 j j+l)(B 1 j-1 j)] 2 = [(AB 1 j)(l j-1 j j+l)] 2 . (33) 

In general, all of these terms contain both physical as well as spurious poles. Physical poles are 
denominator factors of the form (AB i i+1) and (i i+1 j j+l) while spurious poles are all other 
denominator factors. We often refer to physical poles as local poles and to spurious poles as non- 
local. A small explanation for the "non-local" terminology is in order. Consider the 5-particle 
amplitude as an example, where there are three terms in the integrand. These three terms are 

(1234) 2 (AB\ (123)fl(145)) 2 



(AB12)(AB23)(AB34)(ABU) (AB12)(AB23)(AB31)(ABU)(ABA5)(AB51) 
(3451) 



(34) 



(ABM) (ABA5)(AB51) (AB31) ' 

The spurious poles are (AB14) and (AB13). The line defined by Z\ and Z% corresponds to a 
complex point, but what makes (AB13) non-local? The reason is that in field theory 1/(AB13) 
could only come from a loop integration, e.g. it is generated by a local one-loop integral of the form 

~d 4 Z c d 4 Z D ] (OD|(512) H (234)) 



/ 



_vol[GL(2)] 



(35) 



(CDAB)(CD51) (CD12)(CD23) (CD3A) ' 

(This is also the secret origin of the non-local poles in BCFW at tree-level.) 

Back to the 5-particle example, (AB14) and (^41?31) occur each in two of the three terms 
and they cancel in pairs. Indeed upon collecting denominators we find, after repeated uses of the 
Schouten identity, the result for the sum 

(AB12) (2345) (1345) + (AB23) (1345) (1245) + (AB13) (1245) (3245) + (AB45) (1234) (1235) 

( AB12) ( AB23) (AB34) (AB45) (AB51) ' ^ ' 

This is furthermore cyclically-invariant, albeit in a nontrivial way involving Schouten identities. 

Let us also briefly discuss the 6-particle NMHV amplitude at 1-loop. The full integrand has 
16 terms which differs even more sharply from familiar local forms of writing the amplitude. As 
we will review in the next section, the usual box decomposition of 1-loop amplitudes does not 
match the full integrand (only the "parity-even" part of the integrand); even so, there is a natural 
generalization of the basis of integrals that can be used to match the full integrand in a manifestly 
dual conformal invariant form. Any such representation, however, will have the familiar form 
"leading singularity/ Grassmannian residue x loop integral". However, this is not the form we 
encounter with loop-level BCFW. Instead, the supersymmetric 77-variables are entangled with the 
loop integration variables in an interesting way. For instance, one of the terms from the forward 
limit contribution to the 6-particle NMHV amplitude integrand is the following, 

77i(ABl|(23)n(456)) + r/ 2 (4561)(AB31) + 773 (4561) (AB12) 

+ r^ffl (123)0(561)) + V5 (AB1\ (46)0(123)) + 77s (^1 1 (123)0(45)) 
(4561) ( AB45) ( AB61) (ABU) (AB23) (AB13) (ABA1) (AB\ (123)0(456)) (AB\ (123)0(561)) 
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The full expression is given in appendix [Aj Note the presence of the explicit ^independence in the 
argument of the Fermionic 5-function. Seemingly miraculously, when the residues of this integral are 
computed on its leading singularities, the 77-dependence precisely reproduces the standard NMHV 
R- invariants. Of course this miracle is guaranteed by our general arguments about the Yangian- 
invariance of these objects. 



4.4 Unitarity as a Residue Theorem 

The BCFW construction of tree- level amplitudes make Yangian-invariance manifest, but are not 
manifestly cyclic- invariant. The statement of cyclic- invariance is then a remarkable identity between 
rational functions. Of course one could say that the field theory derivation of the recursion relation 
gives a proof of these identities, but this is quite a circuitous argument. One of the initial striking 
features of the Grassmannian picture for tree amplitudes was that these identities were instead 
a direct consequence of the global residue theorem applied to the Grassmannian integral. This 
observation ultimately led to the "particle interpretation" picture for the tree contour, giving a 
completely autonomous and deeper understanding of tree amplitudes, removed from the crutch of 
their field theory origin. 

In complete analogy with BCFW at tree-level, the BCFW construction of the loop integrand is 
not manifestly cyclically-invariant. Again cyclic-invariance is a remarkable identity between rational 
functions, and again this identity can be thought of as a consequence of the field theory derivation 
of the recursion relation. But of course we strongly suspect that there is an extension of the 
"particle interpretation" picture that gives a completely autonomous and deeper understanding of 
loop amplitudes, independent of any field theoretic derivation. 

Just as at tree-level, a first step in this direction is to find a new understanding of the cyclic- 
invariance identities. To whit, we have understood how the cyclic-identity for the 1-loop MHV 
amplitude can be understood as a residue theorem; we very briefly outline the argument here, 



deferring a detailed explanation to 63 . The idea is to identify the terms appearing in the MHV 
1-loop formulas as the residues of a new Grassmannian integral. All the terms in the MHV 1-loop 
formula can actually be thought of as arising from J d 3 ^ZAd 3 ^ 4: ZBY n+ 2,k=2(ZA, Zb, ■ ■ ■), where 
Y n+ 2,k=2 is computed from the G(2,n + 2) Grassmannian integral. Note that Za,Zb appear in 
the delta functions of the integral in the combination C^aZa + C^bZb, so the GL(2)-action on 
(Za,Zb) also acts on (C/3a,C/3b)- Performing the t\a,b and GL(2)-integrals leaves us with a new 
Grassmannian integral: 

.2x(n+2)^ 5 4 (C/3iZj + Cp A ZA + Cp B Z B )(AB) 2 

d C(3a (12)(23)-..(nl) ' (37) 

By construction, this integral has a GL(2)-invariance acting on columns (A,B) and (Za,Zb), and 
hence all of is residues are only a function of the line (ZaZb)- In particular all terms appearing 
in the MHV 1-loop formula, after GL(2) integration, are particular residues of this Grassmannian 
integral. 

As we will discuss at greater length in [63], the equality of cyclically-related BCFW expressions 
of the 1-loop amplitude follows from a residue theorem applied to this integral. In fact, it can be 
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shown that the only combination of these residues that is free of spurious poles is the physical 1-loop 
amplitude. 

At tree level, the cyclic-identity applied to e.g. NMHV amplitudes ensures the absence of spurious 
poles. The same is true at 1-loop level. Since the BCFW formula manifestly guarantees that one 
of the single cuts is correctly reproduced, cyclicity guarantees that all the single cuts are correct. 
Having all correct single cuts, automatically ensures that all higher cuts — and in particular unitarity 
cuts — are correctly reproduced. Unitarity then finds a deeper origin in this residue theorem. 



5 The Loop Integrand in Local Form 

We have seen that the loop integrand produced by BCFW consists of a sum over non-local terms. In 
order to present the results in a more familiar form, and also as a powerful check on the formalism, 
it is interesting to instead re-write the integrand in a manifestly local way (which will of course 
spoil the Yangian-invariance of each term). We will do this for a number of multi-loop examples 
in the next section, but first we must describe a new basis of local loop integrals which differs in 
significant ways from the standard scalar integrals, but which will greatly simplify the results and 
make the physics much more transparent. 

Loop amplitudes are normally written as scalar integrals^] with rational coefficients. Obviously 
this form can not match the full loop integrand, since scalar integrals are even under parity but 
the amplitude is chiral. Let's consider one-loop integrals to begin the discussion. In the usual way 
of discussing the integral reduction procedure, manipulations at the level of the integrand reduces 
integrals down to pentagons 75 . The final reduction to the familiar boxes uses the fact that the 
parity-odd parts of the integrand integrate to zero. 

We are instead interested in the full integrand, however, and since the amplitudes aren't parity 
symmetric, there is no natural division between "parity-odd" and "parity-even". In fact, for the 
purpose of writing recursion relations, it is crucial to know both. Furthermore, the BCFW recursion 
relation guarantees that the loop integrand is dual conformally invariant and thus most usefully 
discussed in momentum-twistor space. We are then led to construct a novel basis of naturally chiral 
integrals written directly in momentum-twistor space, as we now briefly describe. These issues will 
be discussed at much greater length in [63| . 

Let's look at a few quick examples of local integrals written in momentum-twistor space. We 
have encountered the simplest example already; the zero mass integral at 1-loop 



<rZ A a^Z B 
vol[GL(2)] 



(1234) 2 , . 

(38) 



(AB12)(AB23)(AB34)(AB41) ' 



Henceforth, we will drop the integration measure and only write the integrand. The most general 
1-loop integrand is of the form 

(ABYi) . . . (ABY n -4) 



(AB12)(AB23) ■ ■ ■ (ABnl) 



9 Here we abuse terminology and use the term "scalar", which is appropriate at one- loop, to refer to possibly 
tensor integrals at higher-loop order where the tensor structure is the product of "local" factors, i.e., of the form 
{(AB)e i and {(AB) e (AB) k ). 
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where Y-[ , . . . , Y££ 4 are general 4x4 antisymmetric matrices or 'bitwistors'; with 6 independent 
components. Momentum-twistors make integral reduction trivial. Suppose there are 6 or more 
local propagator factors including (ABj\ji+l) ■ ■ ■ (ABj§j§+l) in the denominator. We can always 

expand all the Y IJl s in a basis of the 6 bitwistors Z^Z^^, . . . , Z^Z^^. Inserting this expansion 
into the integrand, each term knocks-out a propagator from the denominator. Thus we can reduce 
any integral down to pentagons. 

These will contain 5 " AB" factors in the denominator and a single " AB" factor in the numerator. 
In the literature, x-space loop integrals are written with numerator factors like (x — Xj) 2 , which in 
momentum-twistor space correspond to {AB j However, we will find more general numerators 

to be more natural. For instance, a typical pentagon integrand we consider takes the form 

(ABU) (5123) (2345) 



(AB12)(AB23)(AB34)(AB45)(AB51> ' 



We can trivially translate this integral into x-space; the numerator is proportional to (x — X14) 2 , 
where X14 is a complex point associated with the line (14) in momentum-twistor space; specifically, 



the pentagon- integral (40) is given by 



< 14 >< 23 > fc 7 „, Wlth XM -\^)- i^n) 



(12)(34)7 (x-xi) 2 (x-x 2 ) 2 (x-X3) 2 (x-X4) 2 (x-x 5 ) 2 ' (14) 

(41) 

The complex point X14 is null-separated from xi, X2, X4 and X5; the second point sharing this property 
is its parity conjugate which will be described shortly. These complex points play an important 
role in the story, and it is most convenient to discuss them on an equal footing with the rest of the 
points by working directly with momentum-twistor space integrands. 

Notably, unlike standard scalar integrals, this pentagon integral is chiral. Like any pentagon 
integral, it has 5 quadruple cuts and twice as many leading singularities. But unlike a generic 
pentagon integral, with this special numerator, half of the leading singularities vanish, and the 
others are all equal up to sign — hence, we say that this integral has "unit leading singularities" . All 
of the local integrals we consider have this quite remarkable feature. 

Local momentum-twistor space integrals can be drawn in exactly the same way as familiar 
planar integrals in x-space; we introduce a new bit of notation to denote the numerator factors. 
The pentagon integral we just discussed is drawn as, 



(42) 



where the dashed line connecting (1,4) denotes the numerator factor (ABM). We will also have 
recourse to use the parity conjugates of these lines. The point Zi in momentum twistor space is 
naturally paired with its projective-dual plane W, = i and the parity conjugate of a line 
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{ij) is the line which is the intersection of the corresponding planes (ij) = i i+1) f] (j— 1 j j+1). 
The numerator factor (ABij) = (A i—1 i i+l)(B j—1 j j+1) — (B i—1 i i+1) (A j—1 j j+1) will be 
denoted by a wavy-line connecting i,j. 

With this notation we can nicely write the integrand for n-particle 1-loop MHV amplitudes as 



/ 



1 

n 



n 



+ 




\ 



+ cyclic 



(nl23)(12 i i+1) 
\ 2 < i < n 



^ , i-1 

(2 j i-1 i) 

x(AB\(123)f)(j-lj j+1)) 
3 < i < j < n 



(43) 
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In this expression we sum over all cyclic integrands, including duplicates, which is related to the 
presence of the 1/n pre- factor. 

For definiteness, we have indicated the numerator factor beneath the corresponding picture. 
Recall the familiar form of the MHV amplitude as a sum over all 2-mass easy boxes; it is amusing 
that in our form the only boxes are 2-mass hard. The algorithm by which this form was deduced 
will be explained shortly. 

We pause to point out that the full integrand for some MHV amplitudes have been computed in 
the literature, in the context of using the leading singularity method to determine the integrand |57] . 
A peculiarity in these papers was that the set of integrals that were used to match all the leading 
singularities did not appear to be manifestly dual conformal invariant — which is particularly ironic, 
given that the leading singularities themselves are fully Yangian- invariant! This led some authors to 
the conclusion that the parity-odd parts of the amplitude are somehow irrelevant, since they not only 
integrate to zero on the real contour but are also not dual conformal invariant. Of course, nothing 
could be further from the truth: we have seen very clearly that the full integrand is determined 
recursively and exhibits the Yangian symmetry of the theory; the decomposition into parity even 
and odd parts is artificial. The problem is quite simple, the basis of scalar integrals has only parity 
even elements! Therefore, one is trying to model the full integrand with a very inappropriate basis. 

From the momentum-twistor viewpoint, the source of the previous difficulties can be seen quite 
explicitly. We have seen that all 1-loop integrals can be reduced to pentagons, but these are tensor 
pentagons, i.e. with factors of AB in the numerator. Now, it is possible to further reduce a pentagon 
with numerator (ABY), with Y corresponding to a real line or not, to a scalar pentagon integral, by 
expanding Y in a basis of the 5 bitwistors appearing in the denominators, together with the infinity 
twistor I IJ . But this breaks manifest dual conformal invariance! Thus the integrands obtained 



in |55,57,59 are indeed dual conformal invariant, but the symmetry was obscured by insistence to 



use scalar integrals. 
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Let's give an example of an interesting two-loop integrand using our notation: 

(1345) (5613) (AB46)(CD\ (234) f] (612)) 
(OD61) (CT>12) (CL>23) (CD3A) (ABCD) (AB34) ( AB45) (AB56) (AB61) 

which we draw as 



(44) 




At two- loops, there are generally 4 solutions to cutting any eight propagators, and so this integral 
has 9 x 4 = 36 different (non-composite) leading singularities. However, the integral is maximally 
chiral: putting any choice of eight propagators on shell will have only one solution with a non- 
vanishing residue. Moreover, the non-vanishing residues are equal up to a sign. This non-trivial 
fact can be understood as following from the global residue theorem applied to the integral. All the 
tensor integrals we write in this paper are chiral in this sense, and the overall normalization of each 
has been chosen so that all its non- vanishing leading singularities are equal to ±1. 

These chiral momentum-twistor integrals have another remarkable feature: they are less IR- 
divergent than generic loop integrals; indeed, many of them are completely IR-finite. Infrared 
divergences arise when the loop momenta become collinear with the external momenta pj . In the 
dual co-ordinate space, this happens when a loop-integration variable x lies on the line connecting 
Xj and Xj + \. In momentum-twistor space, this corresponds to configurations where the associated 
line (AB) passes through the point Zj while lying in the plane (j— 1 j j+1). An integral is IR-finite 
if the numerator factors have a zero in the dangerous configurations. There are an infinite class 
of IR-finite integrals at any loop order; for instance, it is easy to see that the two-loop example 
above is IR-finite. Further discussion of these objects and their role in determining IR-finite parts 
of amplitudes like the remainder [19] and ratio 76 functions will be carried out in [63] . Of course 



we expect that IR finite quantities, such as the ratio function, are manifestly finite already at the 
level of the integrand. 

It is interesting that the naively "hardest" multi-loop integrands can be reduced to finite integrals 
plus simpler integrals. Consider for instance a general double pentagon integrand for six particles, 
of the form 

(ABY 1 )(CDY 2 ) 

(CD6l)(CD12}(CD23)(CD3A}(ABCD}(AB34:}(ABA5}(AB56}(AB61}' 1 ' 

We can expand Y\ in terms of the 6 bitwistors (Z3Z4), (Z4Z5), (Z^Zq), (ZqZi) as well as the 
bitwistors corresponding to (46) and its parity conjugate (46). Similarly we can expand Y2 in 
terms of (Z1Z2), (Z2Z3), (Z3Z4), (ZqZ\) as well as (31) and (31). Doing this reduces the integral to 
finite double-pentagon integrals, plus simpler pentagon-box and double-box integrals. 

Finally, let us describe the general algorithm which we used to find local forms of the loop 
integrands. The first step is to construct an algebraic basis of dual conformal-invariant integrals, 
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over which the integrand is to be expanded. It turns out, quite remarkably, that for at least 1- and 
2-loops an (over-complete) algebraic basis can be constructed which contains exclusively integrals 
with unit leading singularities, in the sense just defined. We have explicitly constructed such a bases 
at 1- and 2-loops and arbitrary n [63] . The second step is to match the integrand as generated by 
equation (29) with a linear combination of the basis integrals. Since the loop integrand is a well- 
defined function of external momenta and loop momenta, this can be done by simply evaluating 
it at sufficiently many random points. Numerical evaluation of the integrand is itself quite fast. 
Finally, this procedure is greatly facilitated by the fact that, when using our particular integral 
basis, the coefficients are guaranteed to be pure numbers (or multiple of leading singularities, for 
arbitrary N fc MHV), as opposed to arbitrary rational functions of the external momenta. 



6 Multi-Loop Examples 

The recursion relation for loops gives a completely systematic way of determining the integrand for 
amplitudes with any (n, k, £). All the required operations are completely algebraic and can be easily 
automated. In this section we use the recursion relation to present a number of multi-loop results. 

As we have stressed repeatedly, the individual terms in the BCFW expansion of the loop 
integrand have spurious poles and are also not manifestly cyclically- invariant; thus as a very 
strong consistency check on our results, necessary for a local form to exist, we verify that the 
integrand is free of all spurious poles: the only poles in the integrand should be of the form 
(i—l i j—1 j), {{AB)i j—1 j), ((AB)£ 1 (AB)£ 2 ) . We also explicitly check cyclic-invariance. Recall 
that the absence of spurious poles and cyclicity guarantees that all single-cuts of the amplitude are 
reproduced, and thus all cuts are automatically correctly matched. While preparing this paper we 
have explicitly checked that our recursive determination of the integrand passes these checks up to 
14 pt N 4 MHV amplitudes at 1-loop, 22-pt MHV amplitudes at 2-loops, 8 pt NMHV amplitudes at 
2-loops and 5-pt MHV amplitude at 3 loops. 

We can expand the integral in a local basis of chiral momentum-space integrals with unit leading 
singularities using the algorithm briefly described in the previous section. While the BCFW form 
of the integrand is almost always more concise than the local form, the local form is more familiar, 
so we will present the results in this way. Indeed, the (modestly) non-trivial work here is only in 
determining the natural basis for local integrands. While this is a straightforward exercise using 
momentum-twistor machinery, the result is non-trivial, yielding a canonical basis of multi-loop 
integrals, which we have constructed explicitly for all n up to 2-loops. In order to present a tree- 
loop result, we also found the 5pt basis at three- loops, deferring a complete discussion to (63] . Given 
the basis of local integrals with unit leading singularities, generating the integrand and finding its 
expansion in the basis is not difficult. The natural basis is over-complete and so the results can be 
expressed in a number of equivalent forms. We will choose the forms that seem canonical and reveal 
patterns. As we will see, somewhat surprisingly, the local forms are also often remarkably simple. 
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6.1 All 2-loop MHV Amplitudes 

The two-loop amplitude for 4- and 5-particles is given by, respectively, 



and 



+ cyclic 
(no repeat) 



3 2 

(2341)(3412)(4123) 



(46) 



+ 




"A 4 

3 2 

(2345)(5123)(3412) (3451)(4513) 

x(AB|(512)n(234)) 



2 + cyclic 

(no repeat) 



while the 6-particle amplitude is 



\ V ✓ 



3 2 

(2345) (6123) (3412) 




+ 



(3456) (4563) 
x(AB|(561)n(234)) 

6 



+ 




1 + 



(47) 




(2345) (3462) 



(3456) (4562) 



+ 



+ 4 



x (AB\ (561) fl (123)) x (AB\ (561) fl (123)) 



+ cyclic 
(no repeat) 




(48) 



(3456)(6123)(4512) 



3 2 

(6235) 
x(AB| (234)0(456)) 
x(CD\ (561) (1(123)) 



To be completely explicit, we have written the numerator factors accompanying each given term 
under its corresponding picture. 

What about higher-points? The parity-even part of the integrand has been computed in |54| , 
though the expressions are lengthy and do not expose a discernable pattern. However, looking 
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at the full (non-parity invariant) integrand for 4-, 5- and 6-particles in momentum-twistor space 
reveals a clear pattern: the structure looks like the "square" of the 1-loop objects, with double-box, 
pentagon-box and double-pentagon topologies. This motivates a simple conjecture for all 2-loop 
MHV amplitudes: 



i+l 



+ 



(n 1 2 3)x 
(1 2 i i i+l i+2) 

2 < i < n 




+ 




fc+1 



(2 j i-1 i)(i-2 i-1 i i+l) 
x(AB\(123)f)(j-lj 3+1)) 

3 < i < j < n 



(2 i j k) 
x(AB| (123)0(^-1 k k+1)) 

x{CD\(i-iii+i)nu-ij j+i)) 

2<i<j-l<k-l<n 



We checked numerically that this matches the 2-loop MHV integrand as calculated by BCFW 
directly. Because the recursion relations are easily automated, this can be verified for any number 
of particles. We have checked this explicitly for up to 22 particles. It is worth emphasizing that 
independent of verifying the local- ansatz, the cancellation of spurious poles (and propagators) is 
a particularly strong consistency check for the recursion relations. For instance, for the 22-point 
2-loop MHV amplitude, there are exactly 49, 590 terms in the BCFW recursion, each riddled with 
spurious poles that cancel in the sum. Even a single sign-mistake would have spoiled this miracle. 

It is interesting to note that the naively "hardest" integrals that appear here — the double 
pentagons — have a numerator which renders them completely finite. 



6.2 2-loop NMHV Amplitudes 

Although structurally identical to the 2-loop 5-particle MHV amplitude, it is worth writing explicitly 
the 2-loop 5-particle NMHV amplitude; it is, 



[12 34 5] 




\ 



V 



3 2 3 

(2345)(5123)(3412) (2345)(3451)(4512) 

x(AB31) 



2 + cyclic 
(no repeat) 



(50) 



Notice how this answer highlights the role played by parity: equations (50) and (47) differ only by 



/ 



the parity of the numerator in the tensor-integral — and one can be obtained from the other simply 
by exchanging wavy- for dashed-lines. 
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Table 1: Coefficients of residue (1) = [2 345 6]. Here, "<?" rotates each figure by g : i i— > i+1, and P 
exchanges wavy- and dashed- lines (together with each figure's corresponding normalization). 




(6234) (6245) 
x(AB53>(CD|(123)n(561)> 



1 + 9 3 

-a(i- fl )(i-P) 




(l + g 3 P) 





(4561) 
x(AB| (345)0(561)) 
x(AB| (612)0(234)) 



(3456) 3 
x (AB| (123)0(345)) 
x(AB\ (456)0(612)) 



(l + g 3 P) 




-(l + g 3 P) 



(3456) (4563) 
x{AB\ (561)0(234)) 




(l+g 3 P) 
x(l+9-9 3 ) 



(2345) (3462) 
x(AB|(561)Q(123)) 




(3456) (4562) 
x(AB|(561)Q(123)) 



(1+9 



1 1 

2 +9 4 ) Y / |(l+9 2 +S 4 ) \ <» / 



(3456) 2 (4512)* 



\ 



3 (2345)(3412)(6123) 2 



(3456) (4512) (6123) 



Next we present the 6-particle 2-loop NMHV amplitude, written in the manifestly-cyclic form, 

(1)/! + cyclic, (51) 

where (1) is the Grassmannian residue given by the R- invariant [2345 6] written explicitly in 
equation ([6]). Below, we show the coefficient I± of residue (1). 

We next move to the 7-particle NMHV amplitude, which will be presented in the form, 



[(7)(l)/ 7 ,i + cyclic] + [(7)(2)J 7j2 + cyclic] + [(7)(3)/ 7 , 3 + cyclic] 



(52) 



where is the Grassmannian residue given by the i?-invariant defined by the complement of 

{i, j} in {1, 2, . . . , 7}. The expressions for J 7) i, I 7 ,2,-f 7 ,3 are given in appendix [B] 
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6.3 3-loop MHV Amplitudes 



r) 



The four-point three-loop amplitude is given by the cyclic-sum of the following two classes of 
integrands: 



+ 




+ cyclic 

(no repeat) (53) 



(2341) 3 (3412) 



(2341)(3412) 
x(AB\ (412) fl(123)) 

Although perhaps visually unfamiliar, the second integral above is commonly referred to as the 
"tennis-court" because of the way it is usually drawn. We have drawn it the way we have to 
highlight the presence of the pentagon sub-integral and the role played by the tensor-integral's 
numerator (which should be read as connecting to vertices "1" and "2"). 

Finally, we give the integrals contributing to the full 3-loop MHV amplitude for 5 particles. It 
is given by the following cyclic-sum of^he integrands, 
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(3451) 3 
x(AB|(234)fl(512)) 



(5123) (4512) (3451) 
x(AB| (123)0(345)) 
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x(AB|(345)f|(123)) x(AB|(451) (512)) 

x(AB| (345)f|(123)) 




(5123) (4512) (3451) 



V 



(5123) (3451) (2345) 
x<^B|(123) 0(451)) 



+ 2 



(3451) (4512) (1234)/(5123) (2345) (3451) /(4512) 




:{AB\ (345) 0(512)) 



x(AB\ (123) 0(451)) 
x(CT>|(234) 0(512)) / 



here, r is the reflection operation that maps i i-> (6 — i). Notice that deriving this three- loop 
amplitude using the loop-level recursion requires both the 1-loop 9-particle N 2 MHV integrand, and 
the 2-loop 7-particle NMHV integrand; and so the success of getting a manifestly-cyclic and spurious- 
polc-frcc, local object is an indirect check of the validity of the whole structure at lower-loops and 
higher points. 

We conclude this quick tour of some simple multi-loop integrands by stressing again a remark- 
able feature of all these results. The integrals that appear are special objects with unit leading 
singularities — they are thus the most natural basis of local integrals with which to match the sin- 
gularities of the theory. As a consequence the coefficients are also simple objects: "±1" for MHV 
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amplitudes, and Grassmannian residues with integer coefficients for more general amplitudes. These 
objects should be thought of as the correct building blocks for the local integrand, just as the BCFW 
terms provide the building blocks for the integrand in Yangian-invariant form. As we will discuss 
below, it is also likely that carrying out the integration will yield "simple" results for these classes 
of integrals. 

7 Outlook 

The loop integrand for scattering amplitudes is a well-defined object for any gauge theory in the 
planar limit, and in this paper we have given an explicit recursive prescription for computing it 
to any loop order in N = 4 SYM, in a way which manifests the full Yangian-invariance of the 
theory. This provides a complete definition of perturbative scattering amplitudes in planar N = 4 
SYM, with no reference to the Lagrangian, gauge redundancies or other off-shell notions. Along 
the way, we have also seen a new physical picture for how loops can arise purely from on-shell data, 
associated with removing pairs of particles in a naturally "entangled" way. From this vantage point, 
a number of directions for future work immediately suggest themselves. 

7.1 The Origin of Loops 

A few years ago, the tree-level BCFW recursion relations sat at an interesting cross-roads between 
the usual formulation of field theory, where space-time locality is manifest, and a hoped for dual 
description, where space-time should be emergent. On the one hand, the recursion relations were 
directly derived from field theory — without the field-theoretic motivation, it was hard to imagine 
the motivation for gluing lower-point objects together in the prescribed way. On the other hand, 
the presentation of the amplitude was very different from anything normally seen in field theory. 
The amplitudes could be presented in many different forms, with remarkable identities guaranteeing 
their equivalence. The simplicity of the answers resulted directly from the presence of non-local 
poles. These properties, together with the dual super-conformal invariance of all terms in the BCFW 
expansions, strongly motivated the search for a dual theory which would make these features obvious, 
and which would furthermore give an intrinsic definition of the tree amplitudes on its own turf. 

The Grassmannian duality for leading singularities provides this dual understanding of tree am- 
plitudes in a satisfying way. The Yangian symmetry is manifest (for all leading singularities and 
not just tree amplitudes). The amplitude can be presented in many forms since it is a contour inte- 
gral, with many representatives for a given homology class. The remarkable identities guaranteeing 
cyclic-invariance (together with important physical properties at loop-level) indeed find a new inter- 
pretation as higher-dimensional residue theorems. And finally, giving the contour integral over the 
Grassmannian a "particle interpretation" poses a natural question, intrinsic to the Grassmannian 
picture, whose answer yields the tree amplitude, along the way exposing a (still quite mysterious) 
connection with twistor string theory. We strongly suspect that a generalization of this picture 
exists that extends the duality to only to incorporate loop amplitudes but also explain why loops 
must be computed to begin with. 

Our extension of BCFW to all loop orders puts loop amplitudes in the same position at the 
cross-roads between field theory and a sought-after dual description that tree amplitudes occupied 
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a few years ago. This should set the stage for fully exposing the dual picture, and we have already 
made some inroads to uncovering its structure. For instance we saw that the remarkable identities 
guaranteeing cyclic-invariance of the MHV 1-loop amplitude indeed have an origin as a residue 
theorem in a new Grassmannian integral closely associated to the "master" integral computing 
leading singularities/ Yangian- invariants. The nature of the "seed" for loops, arising from removing 
particles, is also clearly intimately related to the particle interpretation, which has already played 
a central role in the emergence of locality at tree-level. 

Along these lines, here we give another presentation of the 1-loop MHV amplitudes, which differs 
from the form we obtained using the recursion relation. Consider the tree-level N 2 MHV amplitude 
M n jfc = 2 (-2^1 , . . . , Z m Za, Zb)- The 1-loop MHV amplitude arises directly from the entangled removal 
of A and B: 




Here it is easy to see, using the BCFW form of the tree amplitude, that there is a unique GL(2)- 
contour of integration associated with each term. This formula differs term-by-term from the BCFW 
form of this amplitude. We can however recognize all the terms as residues of the same auxiliary 
Grassmannian integral in equation (37), and we have shown that the equivalence to the BCFW 
form follows from a residue theorem. While this formula does not directly generalize for other 
amplitudes, its form is certainly suggestive. 

Progress on all these questions would likely be accelerated by finding an explicit solution to 
the recursion relation for all (n,k,£), generalizing the explicit solution already known for tree- 



amplitudes 45 



As a final comment, our analysis of loops in this paper has been greatly aided by working 
in momentum-twistor space; these variables allow us to recognize loop integrals in their familiar 
momentum-space setting. However, given that all the elements in the recursion relation were de- 
scribed in manifestly Yangian-invariant ways, it must be possible to translate these results into 
ordinary twistor space. It is likely that the twistor-space formulation will be most fundamental, 
amongst other things it could offer a natural understanding of non-planar loop amplitudes as well. 

The results of this paper also give a renewed hope for extracting loop-information from twistor- 
string theory. As we have seen, loop amplitudes can easily hide in plain sight in subtle ways, mas- 
querading as a formal way of representing "1" in terms of IR-divergent integrals in (3, l)-signature! It 



is likely that a deeper understanding of the contours associated with the "Hodges diagrams" 70 , 77 



already for twistor-space tree-amplitudes in (3, l)-signature, will be important to make progress here. 
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7.2 Simplicity of Integrals and IR- Anomalies 

Putting aside these highbrow issues, we are confronted with a much more urgent question: does 
our understanding of the integrand help us to carry out the integrations to obtain the physical 
amplitudes? Are the symmetries of the integrand of any use? 

In fact the manifestly Yangian-invariant way of presenting the integrand does strongly suggests 
that the integrals themselves will be "simple". The "super" part of super-dual conformal invari- 
ance is already an extremely powerful constraint. Consider MHV amplitudes for simplicity. The 
statement of super-dual conformal invariance is 

K q^j M mkv = 0^-^jM M}1V = for alia, (54) 



where we use the fact that the MHV amplitude has no r\ a dependence. Thus, the only super-dual 
conformally invariant amplitude is forced to be a constant! This reflects the well-known fact that 
the only Yangian invariant with k = is the MHV tree amplitude (=1 in momentum-twistor space). 
Now, we have expressed the integrand for the MHV amplitude in a manifestly super-dual conformal 
(indeed Yangian)-invariant way. Consider for instance the 1-loop amplitude, which has the form 

M MH v = JS A Z A d^Z B F{Z A ,Z B ;Z a ), (55) 

with an entangled contour of integration for Z a ,b] we suppress the explicit expression for F. The 
statement of super-dual conformal invariance is perfectly well-defined at the level of the integrand, 
which turns into a total derivative: 



,f £jM Mm = J SVZ B (rfjjj+dfjjj) F. (56) 



After doing the t]a,b an d GL(2)-integrals, we have 



d ^ mv ,f^^(JLa 1 + JLc %] , ,57) 



Q Z J "my j vo1[ q L(2)] dZ J 

where we suppress the explicit forms of G a AB . We see that super-dual conformal-invariance continues 
to be manifest at the level of the Bosonic loop integrand in the dual co-ordinate space, also at all 
loop orders. 

This symmetry therefore guarantees that no matter how complicated the integrand looks, on 
any contour of integration where the integral is completely well-defined, it can only integrate to 
a constant, "1"! The integral is not "1" only because we choose a contour of integration over 
lines (AB) corresponding to real (3, l)-signature points in dual spacetime, and this integral is IR- 
divergent. We see that IR-divergences are not an annoying side- feature of loop amplitudes, they are 
the sole reason these amplitudes are non-trivial; in this Yangian-invariant form, the loop amplitudes 



are telling us "I diverge, therefore I am" 10 This is a powerful statement that should be turned into 
an engine to simplify the computation of the loop integrals. Due to the IR-divergences, the Yangian 
generators will not quite annihilate the loop amplitude, but they should localize the integral to 



We thank Peter Goddard for this remark. 
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the IR-divergent regions of loop momentum-space collinear to the external particles. In the dual 
co-ordinate space, this is the region localized to the edges of the null-polygonal Wilson loop. It 
seems likely that these IR-anomalies fully control the structure of the amplitude. Amongst other 
things, they must lie behind the astonishing simplicity recently uncovered in the structure of the 



remainder function for the 2-loop, 6-particle MHV amplitude 78 . In the same line of thought, 



it is conceivable that there is a very direct link between the Yangian structure we uncovered and 
the very beautiful connections made at strong coupling with integrable systems, Y-systems, TBA 



equations and the Yang- Yang functional 79 , 80 . Already these developments have allowed a bridge 



to weak coupling by computing sub-leading corrections to collinear limits 81-83 



Having said all of this, there is a very important issue that must be addressed to make progress 
in directly computing these Yangian- "invariant" but non-local integrals. The question is of course 
how to handle IR-regularization for these objects. Dimensional regularization has long been the 
preferred method for regulating IR-divergences in gauge theories, but it does particularly vio- 
lent damage to the structure of the integrand, and is not useful for our purposes. Fortunately, 
there is a better regulator, both conceptually and computationally. Physically, the IR-divergences 



are removed by moving out on the Coulomb branch 53 . This gives a beautifully simple way 



to regulate the integrals in momentum-twistor space which is also useful for practical computa- 
tions |84[|85| . With the loop integrand written in local form, one simply deforms the local propaga- 
tors as (AB j—1 j) i — y (AB j—1 j) + m 2 (AB)(j— 1 j). The physics is always four dimensional. The 
ambiguities in this regulator occur at an irrelevant level C(m 2 )(log(m 2 )) p . In particular there are 
no issues with the notorious "//-terms" in dimensional regularization, and we don't encounter the 
ubiquitous e/e effects either. This is clearly the physically correct regularization for our set-up. 

How should we use this regularization to compute the non-local integrals of interest? One can 
glibly regulate all 4-brackets (ABxy) i— > (ABxy) + m 2 (AB)(xy), but this is not physically correct: 
the regularization of the local propagators is reflecting the (local!) masses induced by Higgsing; and 
so it is not clear how the non-local propagators should be regularized. Indeed, we have checked that 
for the 1-loop MHV amplitudes, this very naive regularization of the integrals does not produce the 
standard result. Of course, since the Yangian invariant form of the full amplitude can be expanded 
in terms of local integrals, we can in principle work backwards to see how the correct local regulator 
affects the non-local integrand; the question is whether there is a sensible way of computing these 
non-local integrals directly. We intend to return to these questions in near future. 

We have emphasized that the Yangian-invariant presentation of the loop integrand strongly 
suggests that the integrals should be simple. But as we have seen in a number of examples, even 
the local forms of the integrand, when written in terms of the natural chiral basis of momentum- 
twistor space integrals with unit leading singularities, look surprisingly elegant. In fact, these 
integrals with unit leading singularities should also be "simple" . The reason is precisely that their 
leading singularities are "1" or "0"; these are the only possible values of the integrals on any closed 
contour of integration, independent of the kinematic variables. This means that e.g. d/dZ T a acting 
on these integrals should also be a total derivative with respect to the loop variables, and that 
they too should be localized to regions with collinear singularities. Since these are local integrals 
their regularization is well defined. Indeed, as we pointed out in our multi-loop examples, the 
naively "hardest" integrals are even IR-finite. The integrals for our form of the two-loop 6-pt MHV 
amplitude have been computed analytically for certain cross-ratios by [34], passing all non-trivial 
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checks. The simplicity of these partial results strongly supports the idea that the full amplitude 
computed with these integrals are also simple. 

7.3 Other Planar Theories 

We end by stressing that many of the ideas in this paper are likely to generalize beyond the very 
special case of M = 4 SYM. Since the integrand is well-defined in any planar theory, one can try to 
determine it with recursion relations just as we have done for Af = 4 SYM. In [50], it was argued that 
the single-cuts of the 1-loop amplitude are well-defined for any theory with at least N = 1 SUSY 
(or J\f = 2 in the presence of massive particles), so the BCFW recursion determines amplitudes at 
least up to 1-loop in these theories too, with or without maximal SUSY and Yangian-invariance. In 
non-supersymmetric theories, further progress on these questions will require a better understanding 
of single-cuts. One difficulty is that the naive forward limit of tree amplitudes is ill-defined. It is 
plausible that this is closely related to presence of rational terms in 1-loop amplitudes, which have 
a beautiful and fascinating structure which is strongly suggestive of a deeper origin. 
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A The BCFW-Form of the 1-Loop 6-Particle NMHV Amplitude 

In this appendix we present the BCFW form of the 1-loop 6-particle NMHV amplitude. The result is 
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A note on notation: the expression (AB1| (56)0(234)) refers to (AB1X) where X = (56) fl (234) 
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B The Local 2-Loop 7-Particle NMHV Amplitude 



Here we give the explicit formula for the 2-loop 7-particle NMHV amplitude. We find it most 
convenient to give a formula for M^ny - M trcc M^v- We can expand this in three cyclic classes 
as [(7)(l)C , 7,i] + [(7)(2)C 7 , 2 ] + [(7)(3)C 7)3 ] + cyclic. We give the expression for the coefficients 
CV,1j CV,2, CV,3 in the tables below. Here "<?" refers to the operation i i + 1, and P is a parity 
flip, that exchanges wavy- and dashed-lines (together with their corresponding normalization), and 
r is the reflection operation i i-> (8 — i). 
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Table 3: Coefficients of residue (7)(1) = [2345 6], continued 
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Table 4: Coefficients of residue (7) (2) = [1345 6] 
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Table 5: Coefficients of residue (7) (3) = [1245 6] 
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Table 6: Coefficients of ^4tree (in addition to the 2-loop MHV amplitude) 
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